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ABSTRACT: Goose fatty liver is one of the most delicious and popular foods in the world, but there is no reliable genetic marker for 
the early selection and breeding of geese with good liver-producing potential. In our study, one hundred and twenty-four 78-day-old 
Landes geese bred in Shunda Landes goose breeding faiTn, Jiutai, Jilin, China were selected randomly. The fatty livers were sampled 
each week after overfeeding during a three week period. Polymerase chain reaction-single strand conformation polymorphism and DNA 
sequencing were used to identify single nucleotide polymorphisms (SNPs) of fatty acid synthase (FAS), which is an important enzyme 
involved in the synthesis of fat under both physiological and pathological conditions. Least-squares correlation was established between 
these SNPs and fatty liver weight, abdominal fat weight, and intestinal fat weight of the overfed Landes geese, respectively. The results 
showed that fatty liver weight of geese with EF and FF genotypes (amplified by primer PI) was significantly higher than that of the EE 
genotype (p<0.05), and liver weight of CD and DD genotypes (amplified by primer P2) was significantly higher than that of the CC 
genotype (p<0.05). Different genotype combinations showed different liver weights, and from highest to lowest were ABDD, DDEF, 
DDFF, DDEE, ABEF, ABFF, AADD, and CDEF. Further analysis of DNA sequencing showed that there were two SNPs within the 5' 
promoter region the FAS gene. The geese of EF and FF genotypes carried a change of T to C, and the geese of CD and DD genotypes 
carried a change of A to G The changes of the bases could potentially influence the binding of some transcription factors to this region 
as to regulate FAS gene. To our knowledge, this is the first report of SNPs found within the 5' promoter region of the Landes goose FAS 
gene, and our data will provide an insight for early selection of geese for liver production. (Key Words: Landes Goose, Fatty Liver, 
Fatty Acid Synthase, Overfeeding, Single Nucleotide Polymorphism [SNP]) 



INTRODUCTION 

The association between gene polymorphism and a 
particular trait has been extensively studied in modern 
animal breeding, and single nucleotide polymorphisms 
(SNPs) are one type that is extremely useful for marker 
assisted selection (MAS). 

Landes geese are famous for their good fatty liver 
production, and the average liver weight after overfeeding 
is about 700 to 800 g. We are interested in the fat deposition 
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in the liver of the overfed goose, a complicated process 
involving many genes related to the metabolism of 
carbohydrates, amino acids, and lipids (Zhu et al., 2011). 
Some reports have shown that SNPs in the genes involved 
in the fatty acid synthesis may associate with fat deposition. 
In turkeys, there was a significant association between fatty 
acid synthase (FAS) gene polymorphisms and fatness 
variability (Sourdioux et al., 1999). In chicken, gene 
polymorphisms in fatty acid binding protein strongly 
correlated with abdominal fat deposition (Wang et al., 2001). 
Five polymorphic sites in 5' and 3' regulatory regions of 
chicken myo statin gene were found to be closely correlated 
with bone growth and fat metabolism (Gu et al., 2003). But 
till now, the associations mainly focus on the chicken, there 
is almost no reports of association between gene 
polymorphism and liver fat deposition, especially in the 
fatty liver of the overfed goose. 
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According to the previous studies, FAS plays important ( 
roles in fatty acid synthesis (Shu et al., 2012), and its 

activity is regulated by insulin signaling and glucose (Han ] 
et al., 2009). Consequently, we wanted to explore whether 

there are SNPs in the FAS gene that were associated with < 

fatty Uver weight. In this study, high carbohydrate diet was < 

used to overfeed Landes geese for 3 weeks, at which time i 

the average weight of a fatty liver could reach 800 to 1 ,000 ' 

g, approximately 5 to 10 time than the normal weight. The i 

single-strand conformation polymorphism (SSCP) ' 

technique was used to divide the geese into different ' 

genotypes. Finally, we found two SNPs in the FAS gene that ' 

were significantly associated with fatty Uver weight, one ' 

was in the promoter region, and the other was in the coding ' 
region. The results of our study could provide useful 

information for geese breeding especially for fatty Uver ' 
production. 

] 

MATERIALS AND METHODS 

Experimental animals ^ 

One hundred and twenty-four healthy 78-day-old 

Landes geese were selected for the overfeeding experiment, . 

which was carried out at the Shunda Landes goose breeding j 

farm, Jiutai, Jilin, China. Each goose was labeled with a , 

number on its wing. Before overfeeding, the average weight , 

of the geese was 4.31+0.25 kg (p>0.05). The overfeeding ^ 

period lasted for 21 days. Finally, the geese were ^ 

slaughtered, and the weight of internal and external fat ^ 

deposition in their Uvers was determined. The breeding ^ 
management and experimental procedure of the geese are in 

conformity with the principles provided by the WMA j 
Declaration of Helsinki. 

Contents of measurement j 

The parameters determined included the weight before ] 

overfeeding, weight at one-week and two-week after , 

overfeeding, weight before slaughter, fatty liver weight, ( 

intestinal fat weight, and abdominal fat weight. All of the j 

weights (Mean+standard error) were determined by an ( 
electronic scale with an accuracy of ±5 g. Before slaughter, 

2 mL blood sample was taken from the wing vein of each ; 
goose, and genomic DNA was extracted for further analysis 

Table 1. Primer sequences for amplification of the FAS gene 
Primer pair Amplified region Base composition (5' to 3') 

P^^ 5'UTR 

P2 5'UTR 

Pg Coding sequence 



of FAS gene polymorphism (Liu, 2003). 

Feeding management of the geese 

The pre-overfeeding period was one week, during which 
each goose was fed about 400 g com every day. When the 
overfeeding began, the geese were placed in the cages (four 
geese in each cage). An overfeeding machine was used in 
our study, and operated by the same person each time. The 
geese were overfed four times each day at an interval of six 
hours with 250 to 350 g high carbohydrate diet. The 
temperature in the goose house was about 20°C to 22°C, 
and controlled in as small fluctuation as possible. The 
relative humidity was 50% to 65%, and internal 
environment of the house was kept quiet and well ventilated. 

Experimental feed formula was composed of 97.6% 
com, 0.9% salt, 1% vegetable oil, 0.5% vitamins and 
minerals. This formula contained 7.8% protein and 14.10 
MJ/kg metabolizable energy. The com was boiled for 30min, 
and then the other ingredients were added before mixing 
and stored until use. 

Primer design 

According to the chicken FAS gene sequence published 
in GenBank, nine pairs of primers were designed using 
Primer Premier 5.0 software, and then synthesized by 
Shanghai Sangon Biological Engineering Technology and 
Services Co., Ltd. Finally, the polymorphic variations in the 
polymerase chain reaction (PCR) products were detected by 
the sodium dodecyl sulfate polyacrylamide gel 
electrophoresis technique. The primer sequences of the PCR 
products with SNPs were shown in Table 1 . 

Polymerase chain reaction amplification 

The PCR reaction system was composed of 2 (xL 
lOxPCR buffer, 0.6 nL 10 \imo\lh dNTPs, 0.6 ^L 10 
mmol/L of each primer, 0.5 |iL Taq polymerase, 1.0 |iL 
DNA template, and 14.7 (xL ddH20. The amplification 
conditions were pre-denaturation at 94°C for 4min, then 34 
cycles of denaturation at 94°C for 45 s, annealing at 60°C 
for 45 s, and extension at 72°C for 1 min, and finally 
extension at 72°C for 7 min. 

Single-strand conformation polymorphism analysis 

The 2 (xL PCR products and 6 (xL loading buffer were 



Product length (bp) Annealing temperature (°C) 
394 49 

228 56 

420 52 



F: TGTGGAGAGGTTGAACTGTG 
R: GATAACGTTTCGGTCAGAGG 
F: AGCAGGAATGGGAACAGGAC 
R: AAGTTCGCTCTGTGTGCAGC 
F: TTGACCGCTGTCTTGATA 
R: GACTGGTCGTGTTCTCCT 



FAS, fatty acid synthase. 
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EF EF EF FF EE DD CD CD CD CC 




-1 



AA AE AE AB 



P.i Pi Ps 

Figure 1. Band patterns of the polymorphic polymerase chain reaction products amplified by three pairs of primers. 
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EE geiiotj'pe FF genotype CC genotype DD geiiotv'pe 

Figure 2. Single nucleotide polymorphisms of EE, FF, CC, and DD genotypes, polymerase chain reaction products of primer 1 and 2. 



mixed, and annealed at 99°C for 10 min. Then the mixture 
was cooled in an ice bath for 5 min, and electrophoresed 
through 12% non-denaturing polyacrylamide gel 
(Acrylamide:Bis solution = 29:1; 10 V/cm, 10 h). Finally, 
the DNA bands were detected by silver staining. 

Cloning and sequencing 

The polymorphic products were recovered and purified 
by QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany), 
and then sent to Shanghai Sangon Biological Engineering 
Technology and Services Co., Ltd. for sequencing. 

Statistical analysis 

SPSS version 19 was used for statistical analysis with 
the following least-squares model. 

Yij = ^+Gi+ey 

yij is the fatty liver weight, |^ is the average value, Gi is 
the fixed effect, e^ is the random effect. 

This model is used for the single genotype analysis. 



y.j = n+G,(Pl)xG,(P2)+G,(P8)+ey 

yij is the fatty liver weight, |j, is the average value, Gj is 
the fixed effect for each genotype amplified by the primer 
PI, P2, and P8, Cy is the random effect. 

This model is used for the analysis of the genotype 
combinations. 

Hardy-Weinberg test was also calculated (Supplemental 
data: See e-version for supplement). 



RESULTS 

Different polymorphic sites in different genotypes 

The PCR products amplified by the three of the nine 
primer pairs had polymorphic sites. According to the 
electrophoretic patterns, PCR products amplified by primer 
1 (PI) were determined as three genotypes, EE, EF, and FF. 
For primer 2 (P2), three genotypes were CC, CD, and DD, 
and for primer 8 (P8), there were only two genotypes, AA 
and AB. The electrophoretic patterns of the primers are 
shown in Figure 1, and the sequencing results of different 
genotypes are shown in Figure 2 and 3. 

We used TRSEARCH ver.1.3 software to predict the 
potential regulatory elements and protein-binding sites in 
the 5' promoter region of the goose FAS gene (URL 
http://www.cbrc.jp/research/DSB/tfsearch.html). The difference 
in regulation and transcription factor binding at the 
polymorphic sites before and after the changes of the bases 
are summarized in Table 2. The changes at these bases 
altered the predicted transcription factor binding in the 5' 



C TG AA GTAGCT ACA 



CTGA GTAGCT ACA 




AA genotype AB genotype 

Figure 3. Single nucleotide polymorphisms of AA and 
genotypes, polymerase chain reaction products of primer 8. 
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Table 2. Prediction for the binding site of transcription factors 
with the 5' region of M5 gene 



Position of 


Mutation 


Transcription 


Transcription 


the binding site 


type 


factor 


factor 


Site 1 (-54 bp) 


T^C 


Oct-1 




Site 2 (-545 bp) 


A^G 




MZFl 



FAS, fatty acid synthase; MZF, myeloid zinc finger protein 1 . 



promoter region in the FAS gene. For polymorphic (T/C) 
site at position -54, the regulatory factor Oct-1 was 
predicted to bind when it was T, but when it was changed to 
C, no regulatory factor was predicted to bind. For 
polymorphic (A/G) site at position -545, no regulatory 
factor was predicted to bind when it was A, but when it was 
changed to G, myeloid zinc finger protein 1 (MZFl) binding 
was predicted. 

Least-squares analysis between fatty acid synthase gene 
single nucleotide polymorphisms and internal and 
external liver fat deposition 

The SPSS 19 was used to analyze the relationship 
between FAS gene SNPs and internal and external fat 
deposition in the goose fatty liver. The results are shown in 
Table 3. 

The fatty liver weight of EF and FF genotypes was 
significantly higher than that of the EE-type (p<0.05), so 
the F genotype favored fat deposition. The fatty liver weight 
of CD and DD genotypes was significantly higher than that 
of CC genotype (p<0.05), so the D genotype was favorable. 
The fatty liver weight of AA genotype was higher than that 



of AB genotype, but the difference was not significant 
(p>0.05). 

Relationship between fatty acid synthase polymorphic 
genotypes and fatty liver performance 

Least-squares analysis on 2x2 and 3x3 multiple 
genotype combinations was performed to test the 
relationship between the composite genotypes and fatty 
liver performance. The results are shown in Table 4 and 
Supplementary Table SI (See e-version for supplement). 

From Table 4, among the combinations of two 
genotypes at two sites, the fatty liver weight of the ABDD 
genotype was the highest, up to 950.23 g, while that of the 
AAEE genotype was the lowest, 585.23 g. From highest to 
lowest were ABDD, DDEF, DDFF, DDEE, ABEF, ABFF, 
AADD, CDEF genotypes (all of which exceeded 800 g). 
The abdominal fat weight of DDFF and DDEF genotypes 
was significantly higher than that of AAEE type (p<0.05). 
The intestinal fat weight of CCFF and AAFF genotypes was 
significantly higher than that of DDEE type (p<0.05). 

DISCUSSION 

Single nucleotide polymorphisms in the 5' promoter 
region of the fatty acid synthase gene may influence its 
transcription 

The products of PI, P2, and P8 were aligned with the 
corresponding DNA sequence of FAS in GenBank. A 
change of base T to C was foimd at -54 bp site in the 5' 
promoter region. The sequence of EE genotype was 



Table 3. Relationship between genotypes and the pcrformanee of geese 



Prinici pail" 




P, 






P: 






Ps 


Genotype 


EE 


EF 


FF 


GG 


GD 


DD 


AA 


AB 


Sample number 


47 


30 


47 


44 


14 


66 


39 


85 


Genotype rate 


0.380 


0.240 


0.380 


0.355 


0.113 


0.532 


0.315 


0.685 


Fatty liver weight (g) 


654.21+24.54" 


829.46+34.19' 


775.33+29.56'' 


641.00+25.36" 


780.12+45.43' 


810.31 + 19.35' 


780.25+19.54" 


730.41+10.87" 


Abdominal fat weight (g) 


413.78+13.09" 


418.66+14.15" 


406.59+10.77" 


389.28+9.91" 


406.56+18.13" 


425.81+8.97" 


417.48+12.78" 


404.70+8.34" 


Intestinal fat weight (g) 


346.76+7.93" 


354.04+12.44" 


371.56+10.01" 


351.01+10.03" 


343.57+21.06" 


353.58+8.90" 


350.71+11.53" 


352.31+6.66" 


PGR, polymerase chain reaction; PI, PGR products amplified by primer 1 were determined as three genotypes, EE, EF, and FF; P2, primer 2, three 


genotypes were GC, CD, and DD; P8, primer 8, there were only two genotypes, AA and AB. 








Different letters indicate counterparts in the same row having significant difference between the average value (Mean± standard error) (p<0.05). 



Table 4. Relationship between different genotype combinations and the performance of geese 



Genotype 


Sample 
number 


Fatty Uver 
(g) 


Abdominal 
fat(g) 


Intestinal 
fat(g) 


Genotype 


Sample 
number 


Fatty Uver 
(g) 


Abdominal 
fat(g) 


Intestinal 
fat(g) 


AACC 


16 


615.32±42.53"' 


429.37±15.72'' 


351.87+17.35* 


ABEF 


34 


856.62±34.32«" 


429.12+11.95" 


357.06±12.71" 


AACD 


19 


652.40+41.37"^ 


385.78+ 16.40" 


368.42+13.16" 


ABFF 


20 


820.21+44.72'"^'= 


417.45+18.05" 


357.50+22.34" 


AADD 


3 


813.45+69.41'"" 


430.00+37.84"' 


333.33+43.34"' 


CCEE 


21 


601.15+32.78" 


434.29+15.89" 


341.43+ 10.38" 


AAEE 


16 


585.23+42.53" 


353.13+20.10" 


337.50+13.83"' 


CCEF 


9 


760.33+60.12"*'" 


417.78+19.70" 


322.22+40.54" 


AAEE 


10 


739.05+63.28"'"" 


397.00+30.22"' 


345.00+17.14"' 


CCFF 


10 


710.23+60.12""^" 


419.00+22.38" 


401.00+26.81' 


AAFF 


9 


682.26+60.16"" 


418.88+20.51"' 


386.67+25.11 ' 


CDEF 


29 


809.09+37.14""' 


413.10+14.38" 


367.93+12.51" 


ABCC 


24 


681.32+32.68"'' 


435.16+12.48"' 


360.00+10.70" 


CDFF 


14 


783.32+50.81"'='" 


391.43+20.16" 


340.71+28.98" 


ABCD 


47 


780.43+209.17"'"" 


418.93+10.14"* 


349.36+11.41" 


DDEE 


2 


910.45+127.53"'= 


410.00+40.00" 


300.00+40.00" 


ABDD 


11 


950.23+48.24= 


434.45±25.02"* 


363.64+22.13" 


DDEF 


7 


944.35+45.39' 


460.71+28.19' 


315.71±28.44" 


ABEE 


31 


689.53+34. IS"'' 


410.97+10.34"' 


345.81+8.71" 


DDFF 


4 


915.01465.08"° 


457.50±14.93' 


317.50+31.99" 



Different letters indicate counterparts in the same row having significant difference between the average value (Mean+standard error) (p<0.05). 
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identical with that of wild type, and FF genotype was a 
mutant. A change of base A to G was found at -545 bp site. 
The sequence of CC genotype was identical with that of 
wild type, and DD genotype was a mutant. A single 
nucleotide deletion of base A was found at 3205 in the 
coding region. The sequence of AA genotype was identical 
with that of wild type, and AB genotype was a mutant. 

TATA box and CAAT box were found between -60 bp 
and -28 bp at the 5' end of the goose FAS gene (Kameda, 
1991). They found that this region was very important for 
gene transcription and expression, as the mutation of 
several nucleotides or even one nucleotide led to changes in 
transcription. In our study, two SNPs were found in the 5' 
promoter region. One was T to C mutation at -54 bp, close 
to the transcriptional start site and the TATA box, which 
could function as an important transcriptional element. 
TRSEARCH verl.3 software was used to predict the 
effects of this change on potential transcriptional regulatory 
elements and protein binding sites, and we found that when 
a T was present at -54 bp, Oct-1 binding was predicted, 
while when it was C, no regulatory factor was predicted to 
bind. According to the correlation analysis, this mutation 
had a significant effect on fatty liver weight, consistent with 
Kameda's study. Oct-1 is a member of the pit-oct-unc 
transcription factor family, and expresses in most cells, and 
regulates the expression of a variety of housekeeping genes. 
Lack of Oct-1 may affect the expression and regulation of 
FAS, but further study is needed to uncover the specific 
mechanism. 

The other was A to G mutation at -545 bp. When base A 
was present at this site, no regulatory factors were predicted 
to bind, but when the base was changed to G MZFl binding 
was predicted. The MZF-1 binding sites are present in the 
promoters of several genes. It can achieve transcriptional 
regulation through binding DNA or RNA, or directly 
regulate gene transcription by forming complexes with zinc 
finger proteins. Lack of MZFl activity greatly increased 
cell differentiation in mice (Gaboli et al., 2001). Here we 
found that the presence of the MZFl binding site might 
significantly affect FAS transcription and fatty liver weight. 
According to the strong correlation of the SNPs and the 
fatty liver weight, we concluded that these two sites could 
be used as candidates for MAS of the fatty Uver-producing 
trait. 

Fatty acid synthase gene single nucleotide 
polymorphisms showed correlation with liver fat 
deposition 

A study on commercial turkey breeding showed that a 
restriction enzyme site of Mspl within FAS correlated with 
fat deposition (SourdioiK et al., 1996). Four hundred and 
ten MicrosatelUte markers on 25 chicken chromosomes 
were screened for quantitative trait loci (QTLs), and found 



that there were multiple QTLs within FAS gene (Jennen et 
al., 2004). So, the SNPs we found could locate in QTLs 
which control the fat deposition in the liver. 

Different genotype combination showed synergistic 
effect on fat deposition 

In recent years, research on the correlation of human 
diseases with genetic variation has shown that correlation 
analysis of combining multiple SNPs with the trait is much 
more precise than that of single SNP. Here, the correlation 
analysis between the combined genotypes of multiple FAS 
gene sites and the fatty liver trait showed that the 
combination of two genotypes among three polymorphic 
sites made a significant difference in the fatty hver weight 
(p<0.05). The ABDD genotype was the strongest 
combination of all genotypes, and the average fatty liver 
weight reached 950.23 g, which was significantly different 
from the weight of the other genotypes (p<0.05). The effect 
of the single-site genotype of P8 on the fatty hver weight 
was not significant, but the combination ABDD genotype 
was significantly different from other combination 
genotypes, and this may reflect the synergistic effect of PI 
and P8. From Table 4, we can see that the genotypes with 
higher fatty liver weight were mostly combination of D and 
F alleles, which suggests that these two alleles could 
function synergistically. The correlation analysis between 
the three-site combination genotypes and the fatty liver 
weight revealed that the combination genotypes did not 
show significant differences in fatty liver traits, which may 
due to the small sample number. But the genotypes with 
larger fatty livers all had D and F alleles, further confirming 
that there was synergy between them. The statistical 
analysis of fatty liver weight with different genotype 
combination showed that B, D, and F were the donainant 
alleles, because the fatty Uver weight of their heterozygous 
or homozygous genotypes were significantly or very 
significantly higher than those of other homozygous or 
heterozygous genotypes (p<0.05 orp<0.01). 

The average fatty liver weight of genotypes without any 
of the three advantageous alleles B, D, or F was 
516.70+100.10 g. The average fatty liver weight of 
combination genotypes containing one dominant allele was 
688.50+174.91 g, which was 71.82 g higher than that of 
other combination genotypes containing no advantageous 
alleles. The average fatty liver weight of combination 
genotypes containing two dominant alleles was 
721.80+171.67 g, which was 33.39 g higher than that of 
combination genotypes having one dominant allele. The 
average fatty Uver weight of combination genotypes 
containing three dominant alleles was 827.58+194.64 g, 
which was 105.69 g higher than that of combination 
genotypes containing two advantageous alleles (p<0.01). 
The average fatty Uver weight of combination genotypes 
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having four advantageous alleles was 64.92 g higher than 
that of those containing three advantageous alleles (p<0.05) 
and was extremely significantly (375.80 g) higher than that 
of the combination genotypes containing no advantageous 
alleles. All the above results fully proved that the D, F, and 
B alleles are the dominant alleles for goose fatty liver 
weight, though due to the small sample number. Further 
large scale experiments will be required to confirm the 
synergy between the three alleles. 

The results also showed that the combination genotypes 
not only affected fatty liver weight, but also affected 
abdominal and intestinal fat weight. Among the single 
genotypes, there was no significant difference for each trait, 
but in combination genotypes of two sites, difference in 
each trait between different genotypes was significant 
(p<0.05). It suggested that different sites might exert minor 
cumulative effect. The body fat level at different growth 
periods are positively correlated with liver FAS activity 
(Tian et al., 1996). The results of this study are consistent 
those results. 

In conclusion, our study identified two SNFs in the 5' 
promoter region of Landes goose FAS gene, and correlation 
analysis showed for the first time that these SNFs 
associated with fatty liver weight. The alleles D, F, and B 
were dominant for goose fatty liver weight, and they could 
be used in MAS for the fatty Uver trait. 
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